#ifndef GAUSS_BANDA_H
#define GAUSS_BANDA_H

#include "mat_pq.h"

/*cada fila de la matriz tiene p+q+2 elementos
 * el lugar i,p+q+1 de la matriz corresponde al resultado de la ecuacion i
 * el lugar i,j para j<p+q+1 corresponde al lugar i,j-p de la matriz real
 * O(n*p*q)
 */
vector<double> solve(const mat& m, int p, int q);

#endif
